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0^ ■ ABSTRACT 

'-^ \ Oxygen to iron abundance ratios of metal-poor stars provide information on nucle- 

I ' osynthesis yields from massive stars which end in Type II supernova explosions. Using 

^ . a standard model of chemical evolution of the Galaxy we have reproduced the solar 

to '. neighborhood abundance data and estimated the oxygen and iron yields of genuine SN 

^ ' II origin. The estimated yields are compared with the theoretical yields to derive the 

relation between the lower and upper mass limits in each generation of stars and the 

^ . IMF slope. Independently of this relation, we furthermore derive the relation between 

the lower mass limit and the IMF slope from the stellar mass to light ratio in the solar 

neighborhood. These independent relations unambiguously determine the upper mass 

limit of m„ = 50 ± 10 M© and the IMF slope index of 1.3 - 1.6 above 1 Mq. This 

upper mass limit corresponds to the mass beyond which stars end as black holes without 

ejecting processed matter into the interstellar medium. We also find that the IMF slope 

index below 0.5 Mq cannot be much shallower than 0.8. 

Subject headings: galaxy: abundances — galaxy: evolution — galaxy: solar neighbor- 
hood — stars: supernovae 



1. INTRODUCTION 

The initial mass function of stars (IMF) is a key 
ingredient in modeling the evolution of galaxies. The 
IMF shape together with the lower and upper bounds 
of stellar mass largely influences the resulting abun- 
dance pattern of heavy elements ejected through su- 
pernovae explosions, the total amount of mass con- 
tained in stellar remnants, and the integrated photo- 
metric properties of galaxies (cf. Tinsley 1980). 

As a conventional practice the IMF is assumed to 
be a time-invariant mass spectrum having a power- 
law of the form n{Tn)dTn oc m~^^'^^'dm [mi < m < 
niu), where n(rn)dm is the number of stars in the 
mass interval m to m + dm. The IMF is usually de- 
rived from the observed luminosity function of stars, 
but even the local IMF in the solar neighborhood is 
still controversial ranging between Salpeter's (1955) 
IMF (x=1.35) and Scab's (1986) IMF (x=1.9) for m 
^ 1 Mq. The integrated photoionization rates and 
colors of nearby disk galaxies are fitted by a Salpeter- 
like shallower IMF (Kennicutt 1983; Kcnnicutt, Tam- 
blyn, & Congdon 1994; Sommer-Larsen 1996). Such 
a global IMF for disk galaxies could be compared 
to the solar neighborhood IMF, because theoretical 
arguments indicate that the IMF is originated from 
fragmentation of gas cloud in a way nearly indepen- 
dent of local physics in the gas although the IMF slope 
might change in starbursts (Silk 1977, 1995; Yoshii & 
Saio 1985; Price & Podsiadlowski 1995). 

This controversial situation for determining the 
IMF slope holds also for the Large Magellanic Cloud 
(LMC). Through transformation of the photomet- 
ric stellar data to the theoretical H-R diagram. Hill, 
Madore, & Freedman (1994) have obtained a steep 
IMF {x = 2), whereas significantly shallower slope 
is derived by Massey et al. (1995) using a similar 
method (a;=1.3), by Hunter et al. (1995) from re- 
cent HST observation [x = 1.22), and by Will et al. 
(1995a,b) from the stellar luminosity function (x ~ 
1.1 - 1.2). Such a difference in the reported IMFs for 
either the solar neighborhood or the LMC suggests 
that the previous methods can not enable a precise 
determination of the IMF slope. 

In addition the upper mass limit rriu remains still 
undetermined. The ratio Y/Z {Y: helium, and Z: 
metals more massive than He) is a strongly decreas- 
ing function of progenitor mass. Then, the relative 
change in helium to metal abundance Ay/AZ from 
the big bang to the present time is strongly depen- 



dent on rriu- Maeder (1992, 1993), based on his stellar 
models with the empirical mass loss function, found 
that the observed high ratio AF/AZ in extragalactic 
HII regions (Walker et al. 1991; Pagel et al. 1992) 
can only be matched with ?Ti„ being as low as 20 - 
25M0. This is significantly smaller than the estimate 
of TO„ - 40 - 8OM0 by van den Heuvel (1992) from 
the presence of black hole X-ray binaries and X-ray 
pulsars in high mass X-ray binaries. 

Recent accumulation of observed heavy-element 
abundances for the solar-neighborhood stars offers an 
opportunity of using a new method for the IMF de- 
termination based on the supernova nucleosynthesis 
argument. In a supernova explosion, specific nucle- 
osynthetic features will be revealed, and the features 
averaged over an IMF arc incorporated in the gas from 
which stars of next generations are born. The stars 
below \Mq do not eject the synthesized heavy ele- 
ments into the interstellar medium, but confine the 
heavy elements produced by massive stars when stars 
are formed. From a viewpoint of chemical evolution, 
it does not matter how stars below 1 Mq are dis- 
tributed. The important point is therefore the total 
stellar mass of such low-mass stars. In other words, 
we can assume a single slope IMF over the whole mass 
range and can know the mass fraction below 1 Mq 
from the combination of the lower mass limit m^ and 
the IMF slope x which can reproduce the observed 
chemical properties. 

The mass of ejecta of elements like O, Ne, and Mg 
from Type II supernovae (SNe II) increases strongly 
with progenitor mass. On the other hand, the mass 
of products from explosive burning Si, S, Ar, Ca, and 
Fe remains close to constant as a function of progen- 
itor mass or even slightly decreases in the case of Fe 
(Nomoto & Hashimoto 1988; Thielemann, Nomoto, 
k Hashimoto 1993; Hashimoto et al. 1993a, 1996). 
Thus, the average ratios of 0/Fe, Ne/Fe, and Mg/Fe, 
as they result from a whole population of SNe II in- 
tegrated over the IMF, will also depend strongly on 
vrLu while a weaker dependence is expected for Si/Fe 
through Ca/Fe. 

The 0/Fe abundance ratio of nearby stars is among 
the best observed ones of these ratios (e.g., Barbuy 
1988; Barbuy & Erderlyi-Mendes 1989; Gratton 1991; 
Nissen & Edvardsson 1992; Edvardsson et al. 1993). 
Yoshii, Tsujimoto, & Nomoto (1996, hereafter YTN) 
used the standard chemical evolution model and de- 
termined the O and Fe yields from SNe II. These 
empirical yields are compared with the updated su- 



Table 1 
The chemical evolution models for the solar neighborhood 



Input quantities at Tq ■ 

IT^Tg) \07W9 



ISGyr 

[WHI7 



Calculated metallieity yields 



Pii,o 



PlI.Fc 



Pla.Fc 



1 


0.19 


-0.03 


0.05 




0.25 


0.00 


0.08 


2 


0.19 


-0.03 


0.05 




0.25 


0.00 


0.08 


3 


0.19 


-0.03 


0.05 




0.25 


0.00 


0.08 



8.03x10-'^ 
9.54x10"^ 
8.21x10-3 
9.54x10-3 
8.42x10-3 
9.66x10-3 



3.12x10-"' 
3.71x10-"' 
3.20x10-4 
3.71x10-4 
3.28x10-4 
3.76x10-4 



7.76x10-4 
9.56x10-4 
7.63x10-4 
9.23x10-4 
7.66x10-4 
9.19x10-4 



Note. — fgiTa) is the present mass fraction of the gas. [0/H]g and [Fe/H]g are the present abundances of oxygen and iron in the 
gas, respectively. 

Note. — The parameters used in common are tia = 1.5 Gyr (the lifetime of Type la supernova progenitors), tin = 5 Gyr (the 
timescale of the gas infall), and [0/Fe]n = +0.4 (the oxygen to iron abundance ratio for metal-poor population II stars). 



pernova nucleosyntheis yields (Hashimoto, Nomoto, 
& Shigeyama 1989, 1996; Thielemann, Nomoto, & 
Hashimoto 1996) to derive the relation between the 
lower and upper mass limits {mi, niu) and the IMF 
slope X. [Wang & Silk (1993) derived the mi — x 
relation using a fixed value of ?Ti„ = 60Mq and 
the averaged O yield over those given by Arnett 
(1978), Woosley & Weaver (1986), Maeder (1992), 
and Thielemann, Nomoto, & Hashimoto (1994).] In 
order to determine a unique combination of (m£, to„, 
x), we need another relation among these IMF pa- 
rameters. 

We note that the stellar mass to light M/L ratio 
provides an independent constraint on the IMF. As 
in the case of chemical evolution, we can assume a 
single slope IMF in calculating this ratio, because the 
luminosity from stars below IMq is negligible and 
therefore how stars below IMq are distributed is not 
important. Moreover, as far as an IMF decreases 
towards high masses, the M/L ratio depends sensi- 
tively on X and mi, but only very weakly on rriu if 
m,u SL 30Mq. Thus, we obtain the m,i — x relation 
which reproduces the present value of M/L in the 
solar neighborhood. 

In this paper we construct the chemical evolution 
model which reproduces the major observational fea- 
tures in the solar neighborhood (§2). Thereby we de- 
rive the m,i — X and m„ — x relations from the oxygen 
to iron abundance ratios of metal-poor stars (§3) and 



the mi — X relation from the stellar mass to light ratio 
(§4). These relations are used to determine m,i, m„, 
and X for the local IMF in the solar neighborhood 
(§5). The result of the paper is discussed in §6. 

2. THE CHEMICAL EVOLUTION MODEL 

SNe II produce most of oxygen after a short life- 
time til ^ 10®-'' yrs due to their massive progenitor 
stars, while SNe la produce most of iron delayed by 
a considerably longer lifetime tja ^ 1 Gyrs of their 
less massive progenitor stars. We employ an instan- 
taneous recycling approximation tu ^ t only for SN 
II progenitors, and two nucleosynthesis approxima- 
tions such as (i) yiw/yu,Fe = (Zo/Zfc)o10[o/^°1" = 
constant where [O/Fe]ii=+0.4 for [Fe/H]<-1 (Grat- 
ton 1991), and (u) 2/ia,o/j/iLO ^ 1 (Tsujimoto et al. 
1995). Hereafter we use the effective yield defined 
as Pi = yi/a where yi is the mass of ejected heavy 
element i and a is the mass fraction locked up in stel- 
lar remnants and low mass stars in each generation 
of stars. [We note that this fraction a and a rate 
coefficient i> of star formation always appear as their 
product in the chemical evolution equations.] The use 
of these approximations in the chemical evolution en- 
ables an empirical estimate of pn,o, Pii.Fo, and pia,,Fo 
without knowing the IMF in advance. 

Following the procedure described in YTN, we con- 
struct a simplified model of chemical evolution allow- 
ing for material infall from outside the solar neigh- 



Fig. 1. — Features of chemical quantities in the solar neighborhood. Shown are the models with tja = 1.5Gyr, 
tin = 5Gyr and fg{TG)— 0.19 for fc = 1 — 3. The symbols represent the data taken from various papers. The model 
metallicity distributions and observed one by Wyse & Gilmore (1995) are normalized to coincide with the total 
number of the sample stars used by Edvardsson et al. (1993). The models for fg{TG)~0.25 are indistinguishable 
from those for fgiTo) = 0.19 and are therefore not shown except in the lower right panel. 



borhood. The chemical evolution equations describ- 
ing time variation of gas fraction fg and the abun- 
dances of oxygen Zq and iron Zpe are then solved 
under the boundary condition that these quantities 
at Tg=15 Gyrs should coincide with those observed 
in the present disk. We assume that the star for- 
mation rate is proportional to some power k of gas 
fraction C{t) = '^[fgit)]'' and that the infall rate has 
a modified exponential form with a timescale tin. For 
a possible range of fc = 1 — 3 we fix tja = 1.5 Gyr 
and iin = 5 Gyr, because both the [0/Fe]-[Fe/H] dia- 
gram and the [Fe/H] abundance distribution function 
of long-lived stars are well reproduced by this combi- 
nation of iia and iin (YTN). 

For the present gas fraction, we set fgiTa) = 
0.19 or 0.25 assuming that the surface mass den- 
sity of the gas component near the sun is 10 or 
11.5 M© pc~^, compared with the total mass den- 
sity of 54 or 46 Mq pc~^, respectively (see §4). For 
the present abundances of oxygen and iron, we set 
i[0/B]g, [Fe/H]<,) = (-0.03, +0.05) for fg{TG)=0.l9, 
and (0.00, +0.08) for /g(rG)==0.25, which are ad- 
justed to give the best agreement with the observed 
chemical properties near the sun. 

Using the above inputs in the chemical evolution 
equations, we have derived the values of pn,0: Pu,Fc, 
and pia,Fc- The results are tabulated in Table 1 and 
are shown in Figure 1. The models for fc = 1 — 3 in the 
case of /g(TG)=0.19 are shown in the [0/Fe]-[Fe/H] 
diagram (the upper left panel), the [Fe/H] distribu- 
tion function of long-lived disk stars (the lower left 
panel), the age-metallicity relation (the upper right 
panel) , and the evolutionary behavior of gas fraction 
(the lower right panel). The data taken from various 
papers are shown by symbols. The models shown in 
the lower left panel are normalized to coincide with 
the total number of the sample stars used by Ed- 
vardsson et al. (1993). We note that the models 
for fg{TG)— 0-25 are indistinguishable from those for 
/g(TG)=0.19 in the first three panels and are therefore 
not shown in these panels. In the last panel, however, 
the model for A; = 1 and /g(TG)=0.25 is also shown 
for reference. 

The empirical oxygen yield pii,o and the star for- 
mation rate C{t) derived here will be used to deter- 
mine the IMF in the following sections. 

The empirical oxygen yield pu^o and the star for- 
mation rate C(i) derived here will be used to deter- 
mine the IMF in the following sections. 



3. APPROACH FROM NUCLEOSYNTHE- 
SIS ARGUMENT 

Comparing the empirical O and Fe yields with the 
updated theoretical supernova yields, we can derive 
the lower and upper mass limits {mi, my) as a func- 
tion of the IMF slope index x. These IMF parameters 
are subject to the following two relations: 



PILO 
PlI.Fe 



Mii,o{m,)n{m)dm 



10 Mq 



Mii,pc{m)n{m)dm, 



(1) 



IOMq 



and 



Mii^o{'m)n{m)dm, 



Pu,o 



IOMq 



M^cni(jn)n{m,)dm. 



(2) 



The ejected oxygen and iron masses {Mu^o, -Wii,Fc) 
from SNe II are taken from explosive nucleosynthesis 
calculations for massive stars with m ~ 10 — TOM© by 
Hashimoto et al. (1989) and Thielemann et al. (1990, 
1996). The nucleosynthesis contribution from to = 
8 — IOM0 stars is assumed to be negligible (Nomoto 
1984, 1987; Hashimoto, Iwamoto, & Nomoto 1993b). 
These ejected masses for each main-sequence stars are 
tabulated in Tsujimoto et al. (1995) and are shown in 
Figure 2. It is seen from this figure that the produced 
oxygen increases steeply as the stellar mass increases. 

The calculated oxygen mass is subject to the com- 
bined uncertainties involved in the ^^C(a, 7)^^© rate 
and convectivc overshooting near the end of core He 
burning. The higher rate and smaller overshooting 
result in the smaller C/0 ratio after He burning and 
thus smaller Ne/0 and Mg/0 ratios after C burning. 
The range of combined uncertainties is constrained by 
the requirement that the resultant SN II yields repro- 
duce the solar ratios of Ne/0 and Mg/0. The allowed 
range of oxygen masses thus obtained lies between 
our standard masses and those by Woosley & Weaver 
(1995), where the latter calculation produces ~ 1.3 
times larger oxygen masses than our standard case. 
Taking into account the above uncertainty of the oxy- 
gen product, we calculate two cases for the oxygen 
product, that is, the standard oxygen product ob- 
tained by our nucleosynthesis calculations (solid line) 
and the high-0 case in which we multiply the stan- 



Fig. 2. — Synthesized and ejected oxygen and iron masses from Type II supernovae as a function of the mass 
of their progenitor main-sequence stars which are assumed to have solar abundances. The masses are in units of 
Mq . Considering an uncertainty of the oxygen products we show two cases such as the standard oxygen product 
obtained by our nucleosynthesis calcultaions (solid line) and the high-0 case in which the standard oxygen products 
for each main-sequence mass are multiplied by a factor of 1.3 (dotted line). 



dard oxygen products for each main-sequence mass 
by a factor of 1.3 (dotted line). 

There also exists an uncertainty concerning the 
iron product, but this uncertainty is relatively small 
as far as an IMF decreases towards high masses, be- 
cause (1) the dependence of the ejected iron mass on 
the progenitor mass is very weak and therefore (2) 
the integrated iron product is mainly determined by a 
lower mass range 10 - 20 Mq as constrained by recent 
observations of Type II, Ib/c supernovae (SN 1983N, 
SN 1987A, SN 1990E, SN 1993J, SN 19941; Nomoto, 
Iwamoto, & Suzuki 1995). Therefore we take into 
account only the uncertainty of the oxygen product. 

The remnant masses Mrem are taken as the masses 
of white dwarfs for the initial low-mass stars, and 
as the masses of neutron stars for the initial massive 
stars. The initial-final mass relations for white dwarfs 
and neutron stars are described in YTN. ( For the 70 
Mq star, Mrem = 1-57 Mq, which is not given in 
YTN.) As for the lower mass limit below which stars 
do not experience the mass loss, we set ?7iwd.i=0.7 Mq 
and therefore M^cm='m below ?7^=r7^wd,^■ 

Using the empirical 0/Fe yield ratio in the LHS 



and the nucleosynthesis product masses in the RHS, 
we derive the upper mass limit rriu as a function of the 
IMF slope index x from eq.[l]. The result is shown in 
the upper panel of Figure 3, where solid and dotted 
lines represent the standard and high-0 cases of the 
oxygen product, respectively. For the standard case a 
Salpeter IMF {x = 1.35) corresponds to to„ = 50 Mq 
and a Scalo IMF {x = 1.9) to m„ = 80 Mq, whereas 
these are to„ = 40 and 50M© for the high-0 case. 
We note that an uncertainty of ?7i„ arising from the 
uncertainty in the oxygen products becomes larger for 
a steeper IMF. 

Similarly using the obtained ruu — x relation in 
eq.[2], we derive the lower mass limit mi as a func- 
tion of x. Since we assume a single IMF slope over 
the whole mass range, the lower mass limit derived 
here might not represent the real mi if the IMF slope 
changes below 1 Mq. The derived mg is a parameter 
to determine the total stellar mass for to < 1 Mq, as 
mentioned in §1. The result of the mi — x relation is 
shown in the lower panel of Figure 3 where two thin 
lines correspond to those for the lowest and highest 
values of pu.o in Table 1. Contrary to the to„ — x 



relation, different choices of Mn,o hardly change the 
me — X relation because the IMF steeply decreases to- 
wards high masses so that the RHS of eq.[2] is not 
very sensitive to niu ■ From the obtained m„ — x and 
TUi — X relations, the mass range covers 0.04 Mq - 50 
Mq for a Salpeter IMF and 0.4 M© - 80 Mq for a 
Scab IMF. 

In the present analysis we do not include a possible 
metallicity dependence of the nucleosynthesis yields. 
Oxygen and Fe are primary yields and do not depend 
on seed abundances from previous stellar populations. 
They could only be affected by a metallicity depen- 
dence of stellar structure, which are not expected to 
be large (Woosley & Weaver 1982). More significant 
is the infiuence of metallicity dependent stellar wind 
losses. If the O output from SNe II is caused almost 
entirely by yields from hydrostatic stellar evolution, 
mass loss would cause a drastic change of 0/Fe with 
respect to calculations which neglect this mass loss 
(Woosley, Langer, & Weaver 1993). In particular, 
the oxygen yield integrated over an appropriate IMF 
becomes smaller than that without mass loss, but 
the difference is not very significant (Wang & Silk 
1993; Prantzos, Vangioni-Flam, & Chauveau 1994). 
This effect of mass loss sets only in for progenitor 
masses of 40-60 Mq with metallicities of [Fe/H] > - 
0.4 (Langer 1989; Maeder 1990,1992; Charbonnel et 
al. 1993; Schaerer et al. 1993). Therefore the m„ — x 
relation stays unaltered because our analysis of the 
0/Fe ratio is based on the metal-poor stars before 
such metallicities are attained. On the other hand, 
given the empirical value of Pn^o, the induced reduc- 
tion in the oxygen yield in the numerator of eq.[2] 
requires lowering the stellar remnant mass in the de- 
nominator, which corresponds to an upward shift of 
the me — x relation in Figure 3. 

4. APPROACH FROM THE M/L ARGU- 
MENT 

The total stellar mass M, is very sensitive to the 
lower mass limit me, and the total stellar luminosity 
L* is sensitive to the IMF slope x rather than m,e- 
For the decreasing IMFs with increasing stellar mass, 
the upper mass limit m,u hardly affects either M* or 
L* as far as iriu ^ 30 Mq because such massive stars 
are very few in number and have very short lifetimes. 
This indicates that the (M/L), ratio is essentially 
determined by a combination of me and x. 

Using the rate of star formation C{t) which re- 



produces the solar-neighborhood chemical quantities 
(§2), we calculate the present stellar mass M^(Tg) 
and luminosity Lv*{Tc) according to the equations 
below: 

M.,{Tg)^ f ^dtC{TG-t) 



/•rrit j-m^ 

X < / mn{in)din + / Mrcm{'m-)n{m)dm 

I ./ min(m^ ,mt) •^ nit 



(3) 



and 



Ta 



Lv*{Tg) — / dtC{TG —t) / lv(m)n{m)dm . 

J tm^ J inin{mi^mt) 

(4) 
The turnoff mass rrit corresponding to an age of t is 
given by (Renzini & Buzzoni 1986) 



\og[mt/MQ] = 0.0558[log(i/yrs)]2 



-1.3381og(i/yrs) + 7.764 



(5) 



and tm^ is the lifetime of stars with mass to„. In 
eq. [3] the first and second terms in the braces denote 
the contribution from main-sequence stars and stellar 
remnants, respectively. The stellar mass-luminosity 
relation lv{m,) in the V band, which is used in eq.[4], 
is taken from the theoretical solar-abundance m - My 
relation compiled by Tinsley (1980). This relation is 
consistent with the recent theoretical relation by Tout 
et al. (1996) and agrees with the empirical relation 
by Scalo (1986). In eq.[4] we do not consider a pos- 
sible metallicity effect in the m - My relation nor 
the luminosity from giant stars. This is because the 
majority of long-lived disk stars of our concern have 
metallicities as high as [Fe/H] ^ - 0.1 (see the lower 
left panel of Figure 1) and there is a large uncertainty 
in estimating the luminosity of giant stars. 

The observed {M/Ly)* ratio, which will be com- 
pared with the calculated ratio, is obtained from 
the surface mass density of nearby disk stars and 
their surface luminosity density of only main-sequence 
stars. Since the dynamical mass in the solar neighbor- 
hood is a sum of masses of stars and gas (Kuijken & 
Gilmore 1989,1991), the stellar contribution is evalu- 
ated from more direct measurements of the other two, 
e.g., subtraction of the gas mass from the dynamical 
mass. Various estimates of the dynamical mass from 
analyses of vertical motion of disk stars are converged 
to give Sdyn=46 - 54 Mq pc^^ (Kuijken & Gilmore 
1989,1991; Gould 1990; Flynn & Fuchs 1994). Much 



Fig. 3. — IJppcr and lower mass limits (m„,?n£) in units of Mq as a function of the IMF slope index x. In the 
upper panel shown are the to„ — x relations for the standard (thick line) and high-0 (dotted line) cases of the 
oxygen product. In the lower panel shown are the TO£ — x relations based on the nucleosynthesis argument (thin 
lines) and the mass-to- light ratio argument (thick lines). 



Fig. 4. — Mass to visual light ratio {M/Ly)* for disk stars in unit of {M/Lv)q as a function of the IMF slope 
index x for different oxygen yields (upper thick line for p = 8.03 x 10"'^ and lower thick line for p = 9.66 x 10^'^). 
The observational bounds of (M/Lv)*=3.1-4.5 are shown by the dotted lines. 



larger values of 67 and 



M, 



pc~2 (Bahcah 1984; 



Bahcall, Flynn, & Gould 1992) are attributed to their 
inadequate choice of tracer stars and restrictive anal- 
ysis technique (Kuijken 1995). The mass of gas is 
investigated in detail by Kuijken & Gilmore (1989), 
and their estimate leads to I]gas=13 ± 3 Mq pc~^. 
Straightforward application of this estimate of Sgas 
gives the gas fraction of 19 - 35 %, the upper bound 
of which is unrealistically high. Indeed Young & Scov- 
ille (1991) study the gas fraction for all type of spiral 
galaxies and report that Sb galaxies like our Galaxy 
have a definitive upper bound of 25 %. Accordingly 
we adopt the gas fraction of 19 - 25 % which corre- 
sponds to S*=35 - 44 Mq pc~^ for the surface mass 
density of nearby disk stars. We note that stellar 
remnants, brown dwarfs and giant stars are not sub- 
tracted out of the above estimate of S*. 

Star count analyses in the direction perpendicular 
to the Galactic disk provide an estimate of the surface 
luminosity density of all stars in the V band such as 
Sy=23.3 LvQ pc"2 (Bahcall 1984), 23.8 Lvq pc'^ 
(van der Kruit 1986), and 25.3 Lvq pc~^ (Yoshii, 
Ishida, & Stobie 1987). These values are consis- 
tent with the value of 24 Lvq pc~^, which is eval- 
uated from 2.4/^m. observations assuming V — K = 



3.2 (Fazio, Dame, & Kent 1990). By subtracting a 
contribution from giant stars, these authors have also 
estimated the surface luminosity density of only main- 
sequence stars as 9.7 Lvq pc^^ (Bahcall 1984), 10.2 
Lvq pc^^ (van der Kruit 1986), and 11.4 Lvq pc^^ 
(Yoshii, Ishida, & Stobie 1987), respectively, leading 
to 9.7 - 11.4 Lvq pc~^ for a reliable range. 

We therefore conclude that the observed (M/Lv)* 
ratio in the solar neighborhood is constrained to be 
3.1 - 4.5 {M/Lv)q- Using this empirical ratio and 
theoretical inputs in eqs.[3] and [4], we obtain the 
mi — X relation as shown by thick lines in the lower 
panel of Figure 3. In eq.[3] we tentatively use ?Ti„ = 
50 Mq, but different choices hardly change the result. 
The allowable range of mi — x relation is confined 
between two lines corresponding to (M/Lv)* = 3.1 
and 4.5 {M/Lv)q. 

5. THE IMF DETERMINATION 

We have obtained two mi—x relations which are in- 
dependent of each other; one relation shown by thin 
lines in the lower panel of Figure 3 is from the nu- 
cleosynthesis argument (§3), and another relation by 
thick lines is from the mass-to-light ratio argument 
(§4). It is evident that two allowable regions overlaps 



in a narrow range of x — 1.3 - 1.6. 



Table 2: Derived parameters for the solar neighbor- 
hood IMF 



the best estimate of niu = 50 ± 10 Mq from Figure 3. 

If we adopt the metallicity-dependent yield, the 
mi — X relations, shown by thin lines in the lower 
panel in Figure 3 (§3), move upwards, leading to a 
slightly shallower IMF. 
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Note. — The upper and lower mass limits m„ and nii 
are in units of solar mass Mq. The symbol /bd stands 
for the mass fraction of brown dwarfs below the hydrogen- 
burning minimum mass of 0.08 Mq but more massive than 

me. 

Instead of displaying two independent m,i — x re- 
lations separately, we substitute the nucleosynthesis 
m„ — X and m,i — x relations in eqs. [3] and [4] to calcu- 
late (M/Lv)* as a function of only x. The resulting 
(M/Lv)* — X relations for different choices of empir- 
ical oxygen yield are shown by thin lines in Figure 4. 
Considering the uncertainty of this relation and an 
observed range of {M/Ly)* ~ 3.1 and 4.5 (M/Lv)q 
(dotted lines), we can obtain x = 1.3 — 1.6 which gives 



6. DISCUSSION 

The IMF slope derived here is consistent with the 
slope X = 1.35 - 1.5 inferred from the photoionization 
properties of nearby galaxies (Kennicutt 1983; Kenni- 
cutt et al. 1994; Sommer-Larsen 1996). Much steeper 
IMFs suggested by Scalo (1986) are ruled out because 
they fail to reproduce the 0/Fe ratio of metal-poor 
stars and the (M/Ly)* ratio in the solar neighbor- 
hood simultaneously. 

The upper mass limit of stars corresponds to the 
maximum supernova progenitor mass beyond which 
the stars form black holes without ejecting nucleosyn- 
thesis products into the interstellar medium rather 
than neutron stars, the formation of which is associ- 
ated with the ejection of substantial amounts of metal 
enriched material. While the details of the supernova 
mechanism still await a final explanation, an observa- 
tional constraint from the 0/Fe ratio of metal-poor 
stars allow to put a limit of to„ = 50 ± 10 Mq on 
the progenitor mass where the transition from neu- 
tron star to black hole remnants occurs. This agrees 
with the estimate of ?7i„ = 40 — 8OM0 by van den 
Hcuvel & Habets (1984) and van den Heuvel (1992) 
from the presence of black hole X-ray binaries and 
X-ray pulsars in high mass X-ray binaries. 

The IMF, having a single slope over the whole mass 
range, extends down to mi ^ 0.03 — O.I6M0 which is 
comparable to the hydrogen-burning minimum mass 
(to « O.O8M0). It is very difficult to determine the 
IMF below IM0 because it is hard to constrain the lu- 
minosity function and there are large uncertainties in 
the mass-luminosity relation. The IMF slope might 
be shallower below 1 Mq and 0.5 Mq (Kroupa et 
al. 1993). Taking this possibility into accout we first 
consider that the IMF slope x below IM0 could be 
shallower than that for the massive part of the IMF, 
provided that the mass fraction below IM0 is kept 
unchanged. Since me for x ^ 1 becomes smaller than 
the theoretical estimate of the minimum Jeans mass 
0.007-0.01 Mq (Low & Lynden-BeU 1976; Silk 1977), 
we restrict ourselves to x = 1.1 — 1.2 and estimate 
how much mass is contained in nonradiating objects 
(brown dwarfs). For reference the mass fraction /bd 
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of such objects below O.OSMq but more massive than 
mi is tabulated in Table 2. Secondly we consider an 
IMF that is shallower below m = 0.5 Mq, which is 
suggested from the starcount analysis by Kroupa et 
al. (1993). The result for the IMF slope x = 0.85 
- 1.0 below m = 0.5 Mq is also tabulated in Table 
2. It is suggested from this table that brown dwarfs 
may exist in the solar neighborhood (see Carr 1994 
for a recent review), but we can not estimate their 
exact mass fraction unless the IMF slope below 1 Mq 
is known. 

An independent estimate of the mass fraction of 
brown dwarfs is possible knowing that the dynamical 
mass is the sum of the gas mass, luminous stars, and 
brown dwarfs. We here use Sdyn = 50 ± 4 Mq pc~^ 
(Kuijken & Gilmore 1989,1991; Gould 1990; Flynn 
& Fuchs 1994), Egas = 13 ± 3 Mq pc'^ (Kuijken & 
Gilmore 1989), and i;i"'»=26 ± 2 Mq pc'^ (Bahcall 
& Soncira 1980; Yoshii ct al. 1987; Mera, Chabrier, 
& Baraffe 1996), where S|f™ stands for the stellar 
mass excluding the contribution from brown dwarfs. 
We obtain the mass fraction of brown dwarfs as 30 
±15%, which is comparable to the other estimates for 
re = 1.1 - 1.2 (m < 1 Mq) or a; = 0.85 (m < 0.5 Mq) 
except for the IMF with x=1.6 for the massive part 
in Table 2. This indicates that the IMF slope below 
1 Mq is not much shallower than the slope above 1 

Mq. 

Our conclusion derived here is supported by the 
recent determination of the IMF slope x = 1 ± 0.5 
for very low-mass stars {m ^ 0.6 Mq) by Mera et 
al. (1996) using both the updated mass-luminosity 
relation by Chabrier, Baraffe, & Plez (1996) and 
the luminosity function where a contamination of 
unresolved binaries is carefully taken into account 
(Kroupa 1995). 
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Science, and Culture in Japan. F.M. acknowledges 
the financial support from the Japan Society for the 
Promotion of Science during her stay at the Univer- 
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this paper. 
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